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Abstract 

We compute the renormalized trajectory of (p'^-th.eoicy by perturbation the- 
ory in a running couphng. We use an exact infinitesimal renormalization 
group. The expansion is put into a form which is manifestly independent of 
the scale parameter. 



1 Introduction 



In Wilson's renormalization group |[W71|, |WK74|| , renormalized theories come as 



renormalized trajectories of effective actions. The renormahzation group leaves in- 
variant a renormalized action up to a flow of renormalization parameters. This 
transformation law defines a renormalized theory. It opens a way to a formulate a 
renormalized theory without the detour to a limit procedure. 

We study the (/("^-theory in four dimensional Euclidean space time. A dimension 
parameter D will be kept in the equations to display the dimension dependence of 
scale factors. We use a renormalization group transformation TZ^ which depends on a 
scale parameter L. The scale parameter is equal to the ratio of an ultraviolet and an 
infrared cutoff. We use an exact renormalization group differential equation for the 
L-dependence of effective actions. On the renormalized 0^-trajectory the effective 
interaction depends on L through an L-dependent running coupling g{L) only. A 
natural choice for g is the effective 0^-coupling. On the renormalized trajectory, the 
normal ordered effective potential V{(f),g) (see (|53|)) satisfies the renormalization 
group differential equation 

- (v<p, ^) I V(0, g) = - (V(0, g), V(0, g)) . (1) 

It is the main dynamical equation in our approach. Its most important property 
is independence of L. The renormalization group leaves invariant the renormalized 
trajectory up to a flow of g. This flow is encoded in a /3-function through an ordinary 
differential equation 

L^9{L)=mL)). (2) 

The operator (T>(j), generates a scale transformation. The right hand side of (|I]) 
is a bilinear renormalization group bracket 5(0)) defined in (|55D. It consists 

of contractions between v4(0) and B{(f)), and is independent of g. The 0^-trajectory 
is the solution to (|l|) with 

V{4>, 9) =9^1 d^x0(x)^ + 9^1 d''xd,<P{x)d,(P{x) + 0{g'). (3) 

It is unique in a space of finite solutions. To be precise, it is unique to all orders of 
perturbation theory in g. Solution is here meant as a formal power series in g with 
polynomial coefficients in 0. The criterion of finiteness is that the polynomials are 
given by smooth kernels on momentum space. Precise definitions are given in the 
bulk of this paper. The first order wave function term in (^ is peculiar to the four 
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dimensional case. The coupling follows when (|T]) is expanded to second order in 
9- 

The program of this paper is an iterative solution of (|1|) and of (|^) in powers of 
g. In its course we do not encounter any divergencies. The result is a perturbation 
theory for the 0^-trajectory which is finite to every order in g. This approach was 
proposed in ||Wi96|| using an iterated transformation with fixed L. In this paper we 
investigate the dependence upon smooth variation of L. It gives a refined formulation 
based on the differential equation (^. We restrict our attention to a one parametric 
family of renormalized potentials. It is distinguished by the property that no other 
vertices are included besides those generated dynamically from the first order 0^- 
vertex. The outcome is then a one- dimensional curve, the 0^-trajectory. We write 
it as a curve of effective potentials whose ultraviolet cutoff is rescaled to unity. To 
make contact with the physical world we would have to supply both a physical scale 
A and a value of g. 

We included elementary details to make this paper self-contained. The experts 
in the field are begged pardon and asked to jump directly to what we call local 
perturbation theory. This paper is organized as follows. In section one we formulate 



the setup, make the connection with [Wi96|, and give a derivation of the infinitesi- 



mal renormalization group. Two comments on the ultraviolet limit in the rescaled 
setup and the relation between Green's functions and the effective interaction are 
included. In section two we translate the renormalization group flow into the nor- 
mal ordered representation. The second section closes with a bound on the bilinear 
renormalization group bracket. In section three we perform first an un-renormalized 
global perturbation expansion as a comparison, and thereafter the renormalized 
cutoff-free local perturbation expansion. Explicit second order calculations are in- 
cluded to demonstrate the methods. In the section four, the renormalization group 
equations are solved. We present a weak bound on the large momentum growth 
of the Euclidean Green's functions. It ensures that the iteration is indeed free of 
divergencies. 



2 Renormalization Group 

The setup of this paper is a momentum space renormalization group for a Euclidean 
scalar field 0. 

2.1 Renormalization Group Transformation 

Consider the following renormalization group transformation TZ^, depending on a 
scale parameter L > 1. Let TZl be composed of a Gaussian fluctuation integral, with 
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covariance and mean ip, and a dilatation Sl of ip. Let the fluctuation covariance 
be defined by 



rL{p) = ^mp)-x{Lp)}, (4) 

where x{p) is a momentum space cutoff function. A convenient choice is the expo- 
nential cutoff 

X{p) = exp (-p^) . (5) 

It will be used in the following. Other choices however work as well, for instance 
Pauli-Villar-regularization. The cutoff function's purpose is to suppress momenta 
outside a momentum slice < < 1- Eq. (^) defines a parameter dependent 
positive operator on the subspace of L2(M.^) consisting of functions /(x) with 
zero mode /(O) = 0. Let d/xrz,(C) be the associated Gaussian measure with mean 
zero on field space. Recall its basic property 

/ d/ir,(C)exp{(C,/)} = exp|i(/,ri/)}, (6) 



and consult for example Glimm and Jaffe ||GJ87|| for further information. 

Let the fluctuation integral of a Boltzmann factor Z{(j)) = exp{—V{(j))} be de- 
fined as the average with respect to (|^), shifted by an external field ijj. A convenient 
notation for this average is 

iZ)r„^ = /d/.r,(C)^(^ + C). (7) 

This fluctuation integral can be derived from a multi-scale decomposition of a free 
massless scalar field. Multi-scale decompositions are reviewed by Gallavotti in ||G85 



See also the recent lectures PG95|| by Benfatto and Gallavotti, and references 
therein. The momentum slice L^^ < |p| < 1 can be thought to label a rescaled 
portion of momentum space degrees of freedom. This portion is integrated out in 
(^. The integration of another portion is prepared for by a dilatation Sl of ip- Let 
the dilatation be given by 



The exponent a = 1 — D/2 is the scaling dimension of a free massless field. Anoma- 
lous rescaling will not be considered below. Non-anomalous rescaling applies at least 
to weak perturbations of a free field. 
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The renormalization group transformation is the composition of (|^ with (^. 
For the Bohzmann factor it reads 



nl^'z{^) = {z)^^^s,^ . (9) 

It is not difficuh find renormahzation group flows where the Boltzmann factor ^(0) 
is a polynomial in the field. In scalar field theory the matter of interest are non- 
polynomial flows of the form Z{(f)) = exp{— ^(0)}, where the potential V{(j)) is 
approximately local and polynomial. In the following it is approximately quartic. 
The renormalization group transformation for the potential reads 

7^P°*F(^) = - log (exp(-\/))r^^^^ . (10) 

The superscripts and will be dropped in the following by neglect of notation. 
The below analysis will be done entirely in terms of the potential. The matter of 
stability bounds on the Boltzmann factor will not be addressed. The method will 
be perturbation theory. It is valid in some vicinity of the trivial fixed point V^*(0) = 
0. The renormalization group transformation for the Boltzmann factor is identical 
with the linearized transformation for the potential at the trivial fixed point. The 
linearized renormalization group is responsible for the leading dynamical behavior at 
weak coupling. A number of elementary properties of (0) are conveniently derived 
from (^), which is why the Boltzmann factor is introduced here at all. 

We restrict our attention to even potentials with V {—</)) = V{(j)). Notice that 
field parity is preserved by ([TI1|). Potentials differing by a field independent constant 
will be identified. We could also impose a normalization condition, for instance 
V{0) = 0. To maintain normalization ([ToD then has to be supplemented with sub- 
traction of TZlV{0). Technically this constant is proportional to the volume, infinite 
in infinite volume. Therefore, (p!0D requires an intermediate volume cutoff to make 
sense. We will wipe this technicality under the carpet and leave (|TU|) as it stands. 
This setup is identical with that in |[Wi96|| up to the scale parameter L, which is 
here variable. 



2.2 Semi-Group Property 

The composition of two renormalization group transformations with scale L is equal 
to one renormalization group transformation with scale L^. Moreover, the renor- 
malization group transformation (p!oD satisfies 

nL,nL, = nL,L„ l,,L2>i, (ii) 

linii^i+TlL = id. (12) 
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The renormalization group therefore defines a representation of the semi-group of 
dilatations of with scale factors L > 1 on the space of effective interactions. The 
proof of this semi-group property is 

7^L,7^i,Z(V^) = Jdfir,^{Ci) Jdfir,^{C2) Z{Sl,{Sl,^ + C,) + C2) 

= J d/^5,,r,,5- (Ci) / d/ir,, (C2) Z {SL,SL,^P + Ci + C2) 

= / d/i5^,ri,5^^+ri, (C) Z {Sl.Sl.iP + C) 

= nL,L,Z{tP). (13) 

It uses a dilatation- and a convolution identity for Gaussian measures. Both follow 
from the functional Fourier transform @. It also uses the elementary properties 

(SljTliiSJ^ + = TliLj) (14) 
SliSl2 = <SliL2 (15) 

of and (^. The property (12) follows from limi^i+ Tl{p) = 0"*", turning (^ into 
a functional ^-measure. 

Let us remark that Tl{p) becomes negative for < L < 1. Gaussian inte- 
gration with negative covariance requires further regularity of functionals on field 
space. Non-invertibility of the renormalization group has been recently emphasized 
by Benfatto and Gallavotti in ||BG95|| . Effective theories tend to require less pa- 



rameters because the number of effective degrees of freedom shrinks in the course 
of a renormalization group flow. Inverse renormalization group transformations can 
however be given a meaning on the renormalized trajectory. It is one dimensional 
from the beginning. Furthermore, an interaction on the renormalized trajectory is 
always the renormalization group image of another one. Restricted to the renor- 
malized trajectory, the renormalization group deflnes a representation of the full 
dilatation group. 

Due to the semi-group property the iteration of renormalization group trans- 
formations with flxed scale is identical with an increase of the scale in a single 
transformation. The renormalization group associates with an initial Boltzmann 
factor Z{(f)) an orbit 

Z(0,L) = 7^z.Z(0), (16) 
parametrized by L > 1. Due to the semi-group property it interpolates the sequence 

7^2z(0) = 7^inZ(0), (17) 

obtained by iteration of TZ^. This interpolation is the motive of the present investi- 
gation, in conjunction with the previous work in | Wi96|| on the discrete case. The 



continuous point of view has the advantage to allow for inflnitesimal renormalization 
group transformations, which can be expected close to the identity. 
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2.3 Infinitesimal Renormalization Group 



The infinitesimal renormalization group was invented by Wilson ||WK74]| . Its im- 
portance as exact dynamical implementation of scale transformations is receiving 
increasing recognition. See also Wegner's review ||We76|| . 

The renormalization group orbit ([l6|) satisfies the remarkable functional differ- 
ential equation 



4 



I f s 



Here T> is the generator of dilatations of the field 



V(f){x) = Si 



(18) 

In real space it is given by 

0(x). 



L=l 



\--- — 



(19) 



Furthermore, C is a rescaled scale derivative of the fluctuation covariance. Its ex- 
plicit expression is 



C = Sr 



(20) 



For a general cutoff function of the form x{p) = Fiv^)-, the covariance (|20|) is di- 
agonal in momentum space with eigenvalues C{p) = —2F'{jP'). In the case of the 
exponential cutoff (|20|) , it is C = 2x- In particular, it is independent of L and 
positive. The proof of ( pJ] ) is 



L^Z(^,L) =L-^j d/irJC) ^(V'l + C) 



d 

= L- 



1 
2 



4- 



dtp 



,5 



L-1 



d/ir,(C) Z{^L + 
5 M 



5^ 



[21] 



It uses the change of covariance formula for Gaussian measures and a dilatation 
identity for functional derivatives. Here ifji stands for Siip. The first operator on the 
right hand side of (pT|) performs an infinitesimal dilatation of the field ip. Without 
this dilatation term we would have a functional heat equation. For the following 
analysis it is important to rescale the field and thus carry along the dilatation term. 
We should however concede that the question of ultraviolet and infrared limit can be 
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approached in a non-rescaled formalism as well. There the concept of a renormalized 
trajectory is hidden in scaling properties of a limit theory. 

It is then straight forward to derive a functional differential equation for the 
potential V{(I),L) = — log(Z(0, L)). It is given by 

Here we have collected the linear terms on the left hand side. Notice that the 
linearization of (|22| ) at zero potential coincides with (]T9|). Variations of this flow 
equation have proved to be useful both in the study of perturbation theory and 
in numerical studies of renormalization group flows. Polchinski |[P84|| for instance 
uses a flow equation without dilatation term in his beautiful proof of perturbative 
renormalizability. 

2.4 Ultraviolet Limit 

Eq. (Q) suggests that we are always performing an infrared limit upon iteration of 
(p!0D. In the rescaled formalism ultraviolet and infrared limit are closely related. We 
include this comment on rescaling to prevent confusion at this point. Consider a 
general massless covariance 




with two sided cutoffs. The ultraviolet limit refers to sending the upper cutoff Auv 
to oo in 

Z'diW = 7^A:H„Aw^tv(^) = / d/^^A:H.A.v (0 ^tv(^ + 0, (24) 

keeping the lower cutoff Air flxed at a renormalization scale. Equivalently, the ratio 
L = Auv /Air is sent to inflnity at flxed lower cutoff Air. Deflne a rescaled bare 
Boltzmann factor by 

Z^»-(0) = Z''-^(5a,>) = Z^»-(0). (25) 

The bare Boltzmann factor is here written in units of the ultraviolet scale Auv; to 
be precise in terms of a rescaled fleld (and rescaled couplings). Analogously write 
the effective Boltzmann factor in units of the infrared scale Air, 

ZliiW = Z^"{Sa,^^) = Z^"W. (26) 
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Then the renormahzation group transformation for the rescaled quantities is pre- 
cisely of the form (|), 

= 7^LZ*"'■"(^^). (27) 

Thus keeping the infrared cutoff fixed, the ultraviolet limit is indeed equivalent with 
an infrared limit for the rescaled system. The proof of (^) is 



(c)^a;v('5a,,-^ + c) 



|d/xrJC) Z'^^\SLi^ + C). 



f28) 



In practice it is convenient to put Air = 1 in physical units. Then the outcome of 
the rescaled renormahzation group is already the desired effective potential. The 
argument can be summarized in the following diagram: 



zbare 



Air 



(0) 
Z^ff{il)) 



In the discrete renormahzation group built upon iteration of (P), rescaling can be 
viewed as a stack of these diagrams on top of each other. 



2.5 Green's Functions 

The effective potential considered here is the generating function of free propaga- 
tor amputated connected Euclidean Green's functions. Let us also include a brief 
mention of this property for the sake of completeness. Let W^^^^{J) be the gener- 
ating function of the connected Green's functions with free propagator i^Air.Auv 
vertices Vl^l^cp). Then 

VaI^W = Wl(l{J) + ^ (j, VA,,„A,,J) , ^P = VA,,„A,J. (29) 
Eq. (^) follows from 

exp{-y;//(^)} = / d/i.,^^,,^^(C) exp{-l^^^-;(^ + C)} 

= exp {-^ (^,<,A,vV5)} / d/^^A:R,Aw(0 exp {-V^^v (C) + (C, <,a,v^) } ' 

^ ^ ^ (30) 
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which uses a shift identity for Gaussian measures. The free propagator amputated 
connected Green's functions are given by the kernels of the development 

^1 (2n)! J ^3^^ 
and the non-amputated connected ones by 

= E 7^ y d^Xi . . . d^X2n J(X) . . . J(X2„) H^Xf/,2n(a:i, • • • , 



(32) 



The non-amputated connected Green's functions are reconstructed from the ampu- 
tated ones through 

J d^Vl d^y2 ^^AiR„Auv(^l - yi) ^^AiR,Auv(^2 " ^2) ^Air^2(J'1' J'2), (33) 

and 

^Af/,2n(^l,--- ,X2n) = 
d^yi . . . d^y2n t^AiR,Auv(^l - Vl) ■ ■ ■ ^^AiR,Auv(^2n " y2n) Vl^j!,^2iyi^ ■ ■ • > l/2n)- 

(34) 

The connection between the rescaled amputated Green's functions and the non- 
rescaled amputated ones on the other hand is 

VXrrW = ^^^"W, Sa.J = ^- (35) 
The non-rescaled amputated Green's functions are therefore explicitely given by 

KZ2ni^l, ^2n) = AlR^^^+^Z') V^Sn (AirXi , • • • , AlRX2n) (36) 

in terms of the rescaled ones. In the ultraviolet limit the infrared scale Air is kept 
fixed and (^) amounts to a finite rescaling. But then the dictionary is complete. To 
obtain non-rescaled non-amputated connected Green's functions from the effective 
potential one has to undo rescaling using ( ^6|) and thereafter undo amputation using 
(33) and (34). It is clear that the infrared scale Air is here an extra piece of 
information which has to be supplied from the outside. 



3 Normal Ordering 

We choose to represent the potential in terms of normal ordered products. The 
payoff of normal ordering is that the linear part of the flow equations simplifies to a 
scale derivative and a dilatation term. The price to pay is a more involved bilinear 
term. 
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3.1 Normal Ordering Operator 

We introduce another covariance which will serve as normal ordering covariance. 
Let V be given by 



1 



(37) 



a massless covariance with unit ultraviolet cutoff but without infrared cutoff. It 
satisfies 



Sr-l(v-rr)S 



(3J 



As is shown for instance in |[Wi96|] , this is the fixed point condition for a flow of 
normal ordering. In dimensions D > 2 the infrared singularity of (|37D is integrable. 
This integrability is sufficient for the below purposes, in particular for an estimate 
on the bilinear renormalization group bracket in eq. (p^. Associated with (|37D is a 
normal ordering operator acting as 



(39) 



on polynomials, more generally power series, of the field. The commutator of the 
normal ordering operator (|39D with the dilatation operator in the functional differ- 
ential equation is computed to 



exp 



-16 6 

— -:r-,V — 



6 



6 



6 



exp 



6 6 



(40) 



The intention with the normal ordering covariance (^) was to obtain this identity. 
Differentiation of (|37D with respect to the scale parameter L supplies us with 



Vv + vV^ 



-S 



d 
'dL 



TrAS 



L-1 



-c. 



(41) 



From it we conclude that the commutator of the functional Laplacian, built from 
the normal ordering covariance, with the dilatation operator is given by 



5 6 



6 



-^7, (Vv + vV'^) 



(42) 
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But this is the infinitesimal version of PD|). Eq. (5D) then follows by integration. 



Since the normal ordering covariance ([37| ) is independent of L, the fiow equation 
( p!9D is equivalent to 

This equivalence suggests a normal ordered representation for the Boltzmann factor. 
We write it in the form 

Z{<P, L) =: L) exp | ^ (^^, t;^^ | L). (44) 

We call -2(0, L) normal ordered Boltzmann factor. Strictly speaking it is the pre- 
image of a normal ordered Boltzmann factor by the normal ordering operator (pOf ). 
Therefore it is not decorated with normal ordering colons. The normal ordered 
Boltzmann factor then satisfies the first order functional differential equation 

Z(0,L)=O. (45) 

A noticeable feature of the normal ordered representation is that the exact renormal- 
ization group equation ( ^5)) exactly performs an infinitesimal scale transformation. 
It is equivalent to 

lAz(0^,L)=O, (46) 

with rescaled field 0l = iSl0- In terms of the potential the theory is not as simple 
as this due to the non-linearity of (^2]). We remark that in perturbation theory we 
have in every order to deal with no more than polynomials in 0. Normal ordering 
thus requires no more than a finite number of extra contractions with the normal 
ordering covariance. They will be shown to be finite. 





3.2 Homogeneous Solutions 

The scaling fields of the trivial fixed point are polynomial solutions to (^5]) . They 
are given by homogeneous kernels. Let us have a brief look at them because they 
will be used to parametrize the potential. A detailed discussion of scaling fields can 
be found for instance in Wegner's review |[We76|] . A polynomial 



C(0, L) = —l d^Xi . . . d^Xn 0(Xi) . . . 0(x„) 0„(Xi, ...,Xn,L) (47) 



w. 
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in the field is a solution to the flow equation (^) iff the kernel satisfies 
( d / D\ " (9 1 

+ - E Xmg^\On{Xi,...,Xn,L) =0. (48) 

Eq. (|48|) is a homogeneous scaling equation. It is apparent that the renormalization 
group flow is a pure scale transformation in this formulation. The general solution 
of (^) is 

OniXi, ...,Xn,L)= L^^^^'/^^OniLx^, • • • , LXn). (49) 

A homogeneous kernel of degree k, 0„(Lxi, . . . , Lx„) = L''(9„(xi, . . . , a;„), yields 
an eigenvector of the dilatation generator with eigenvalue n{l + Z^/2) + k. It is 
called a scaling field. The eigenvalue is called its real space scaling dimension. The 
kernels with negative degree of homogeneity have powerlike fall-off at large distances 
and singularities at small distances. Kernels with positive degree of homogeneity 
are non-local and discarded in field theory. A proper mathematical setting for field 
theoretic kernels is that of symmetric distributions on n copies of real space, given 
by Fourier integrals 

On{Xi, . . . , Xn, L) = 
I (^TtT^ ■ ■ ■ 7^71^ J2P"^^"^ {27t)^6^^^ (J2P"^ dn{pi,...,Pn~l,L). 



m=l / \m=l 



(50) 



We will restrict our attention to translation invariant kernels. The 5-function due to 
conservation of total momentum is then conveniently factorized from the momentum 
space kernels. In the sequel it is understood that the nth momentum of an ra-point 
kernel is p„ = — I]^=iPm- Furthermore, it is understood that the momentum space 
kernels are symmetric functions of the n — 1 momenta. Eq. (50) is a solution to 
©iff 

{4 - + " t.'-i:} ■ ■ ■ ^ (51) 

The general solution to (^Tj) is of course given by 

a(pi, . . . L) = f ^, . . . , (52) 
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The power-counting of a homogeneous momentum space kernel of degree k is there- 
fore D + n{l — D /2) — K. A momentum derivative reduces the power-counting of 
a kernel by one unit. The relevant parts of kernels in field theory come as sums of 
scaling fields with positive momentum space scaling dimension. They are extracted 
by Taylor expansion in momentum space. Here field theoretic momentum space 
kernels will be required to be symmetric, Euclidean invariant, and regular at zero 
momentum. 



3.3 Normal Ordered Potential 



Normal ordering reduces the flow of a Boltzmann factor to a pure scale transfor- 
mation. This is not the case for the potential. Nevertheless it is useful to write 
potentials in normal ordered form at least for the case of weak coupling, where nor- 
mal ordered products are eigenvectors of the linearized renormalization group. Let 
us define 



V^(0,L) =: V(0,L) :,= exp 




V(0,L). 



(53) 



in likeness to (pj). This normal ordered potential then obeys the following non-linear 
functional differential equation 




V(0,L) = -(V(0,L),V(0,L)). 



(54) 



It will be the main dynamical equation of this investigation. The non-linearity 
consists of a bilinear renormalization group bracket 



(V(0,L),V(0,L)) = 



(55) 



Here 0^ and 0^ denote two independent copies of 0. The bilinear term consists 
of contractions between two copies of the potential. Notice that every contraction 
contains one hard line C and any number of soft lines v. Eq. ( |54[ ) has been used 
before in a proof of Symanzik's improvement program | Wi88| . 

As a byproduct we obtain a normal ordered functional differential equation for 
renormalization group fixed points. It is given by 



P0, 



6 



V,(0) = (V,(0),V,(0)). 



(56) 
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Renormalization group fixed points are global solutions to (pB]). Non-perturbative 
tools for their investigation would be a major break through in this theory. Presently 
our toolbox contains only the e-expansion and numerical recipes for truncated sys- 
tems. An investigation of (0) along these lines will be presented elsewhere. 

Eq. ( |56D has a trivial solution V*(0) = 0, the free massless field. In the following 
we will restrict our attention to weak perturbations of this trivial fixed point. It is 
then appropriate to use iterative methods to solve (0). 

The proof of (|5^) is 



V{<P,L) 



L^~(v<j),jj]}V{(j),L) 




(57) 



Undoing the field rescaling 



becomes 



L—V{<Pl,L) 



(5J 



Since C is positive definite, the left hand side of (|58|) is always negative. It follows 
that the scaled potential decreases by value under a scale transformation. Notice 
that the difference between the bilinear term in the non-normal ordered formula- 
tion (|2^ ) and the normal ordered one (|5^ ) is a re-normal ordering operator. In 
perturbation theory it can be seen to create additional loops with normal ordering 
covariance. Moreover, normal ordering covariances appear only in contractions of 
vertices and never at external legs. 
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3.4 Bilinear Renormalization Group Bracket 

In the non-hnear formulation the bihnear term is responsible for inhomogeneous 
terms in scaling equations to be considered below. Let us write it explicitely for 
even monomials 02n{4') in the field 0, given by 

C'2n(0) = TTTTT / '^^^1 • • • d^^2n 4>{Xi) . . . 0(x2n) 02n{Xi, . . . , X2„). (59) 



(2n)! 

Recall that the bilinear term does not depend on the scale L. The bilinear operation 



on two monomials of the form (|59D can be decomposed into 



ra+m— 1 

(O2ra(0),O2ra^(0)) = E A^,,ra^,^ ((^2n * 02ra^)2/ (0) , (60) 

l=\n—m\ 

and is itself a sum of monomials 

{02n * C'2rra)2; (0) = 
J d^Xi . . . d^X2l 0(Xi) . . . (f){x2l) {02n * 02m)2l (^1. • • • > ^2l), (61) 

whose kernels are given by a multiple convolutions 

(C^2ra * 02m)2l (^^l; • • • ! ^2l) = 

y dyi . . . d?/2{n+n^-/) C'(2/l - ?/ra+ra^-^+l) \{ ^(Z/fc " ?/ra+m-«+fe) 

) •^n—m+h Vli • • • ) y n+m—l) 
(^2m{Xn-m+l+l, • • • ; 2;2Z, 2/n+m-i+l! • • • ; 2/2(n+m-/)) + 

((2/)! - 1) permutationsj. (62) 

The kernels are understood to be symmetric under permutations of their entries. 
The multiple convolution involves one hard propagator C and n + m — / — 1 soft 
propagators v, which is at the same time the number of loops. Furthermore, (|60|) 
involves a combinatorial factor 

' ' (n + m — / — 1)!(?2 — m + /)!(m — n + /)! ' 

coming from the number of ways in which the contractions can be made. Eq. (|62|) 
can be interpreted as result of the fusion of two vertices. In the process of fusion 
links are created, consisting of propagators. 
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We present an elementary estimate on this fusion product to get an idea of 
what kind of analytical properties can be expected for the effective potentials. The 
estimate uses an Loo^^-norm in momentum space. The estimate works at this point 
for e > 0, not too large. Later it will be used for e > only. Notice to begin with 
that 

||C'||oo-2e, ||^^||l,-2e < OO, (64) 

for the propagators with exponential cutoff. [| At e = we have for instance ||C||oo = 
2 and \\v\\oo < 2tt^^'^/(D — 2). If the Fourier transformed kernels now satisfy the 
bounds 

||C^2n||oo,e7 ||C^2m||oo,e < OO, (65) 

that is, are finite in the Loo,e-normJ^ then it follows that all the terms (|62|) in the 
decomposition of the bilinear operation have finite Loo,e-norms in momentum space. 
They obey 

\\{02n*02m)2l\\oo,e < ^ \\C\\oo,-2e || ^?|| S"'^' \\02n\\o.,e || || oo,.- (66) 

Therefore, the renormalization group flow preserves the Loo,e-norm of momentum 
space kernels for finite scales. It will be shown that the Loo,e-norm is also preserved 
in the iterative solution of (|l]). The estimate immediately follows from the Fourier 
transform 



{02n ^ 02m)2l (Pl^ ■ ■ ■ ,P2/-l) = 

1 p AD„ AD„ n+m-l-1 

2(20! ' (2^ ■ ■ ■ (27r)^ ^ /J^ ""^^^ 

, Pn-m+h Qli ■ ■ ■ 1 Qn+m-l-l) 

n—m 

((2/)! — 1) permutations I . (67) 



iHere ||C||oo,-2. = sup^^^o {\vip)\e^'\P\} and ||^J||i,-2. = (2^)"^/ d^p|w(p)|e2^IPl denote the 
Loo,-2e- and ii^_2e-norms in momentum space. 

^The Loo,€-norm for the momentum space kernels is defined as ||02n||oo,e = 

SUP(p,,...,p,„)g^,„{|0(pi,... ,P2n)|e-'(|P^I+-+|P="l)} with T'sn = { (pi , . . . , P2„ ) S 

|Pi H + P2n — 0} the hyperplane of total zero momentum. 



X 
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of The ^-functions from translation invariance have again been removed. The 
sums of momenta in the kernels are zero through 

2/ — 1 n—m+l n+m—l—1 

P2l = Pm^ Qn+m-l = - ^ Pk - Yl ^k- (68) 

m=l k=l k=l 

The idea with the parameter e is to use part of the exponential large momentum 
decay of the fluctuation and normal ordering propagators to compensate a possible 
large momentum growth of the kernels. In the initial value problem for ( p^ with 
Loo-bounded inital data this might seem unnecessary. For instance a pure 0^-vertex 
is constant and thus Loo-bounded. The evolution preserves Loo-boundedness for all 
finite scales L. However we cannot expect the solution to be Loo-bounded uniformly 
in L. The limit L — oo requires a separate treatment of zero momentum derivatives 
and Taylor remainders of the non-irrelevant kernels. The price to pay is a growth 
in momentum space. In the four dimensional case it is polynomial in powers and 
logarithms of momenta. With an exponential bound we are very far on the safe side. 



4 Perturbation Theory 

In perturbation theory the effective potential comes in form of a power series 

oo r 

V(0,L,^?) = 5:^V«(0,L) (69) 

r=l ' • 

in a coupling parameter g. The effective potential is here assumed to be zero to 
zeroth order. The effective potential is thus expanded around the trivial fixed point. 
In the sequel (^) will be treated as a formal power series in g. The important 
question of non-summability of ( |69D will not be addressed. 



4.1 Global Perturbation Theory 

We speak of global perturbation theory when the expansion parameter is indepen- 
dent of the scale L. This expansion is not appropriate in the limit L — oo. Divergent 
terms appear and call for renormalization. We nevertheless develop the global ex- 
pansion to some detail to see power counting at work. Inserting (KSf) into the bilinear 



bracket (|55D and organizing the result again in powers of g we obtain 



oo r r—1 / \ 

L, g), V(0, ^, ^7)) = E 7! E □ (V^^H^, L)M-^\<j^, L)) . (70) 



r=2 ' ■ s=l 



Let us introduce the abbreviation }C^^\(j),L) for the sum of brackets to order r 
on the right hand side of (^). The effective potential (^) is a solution to the 
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renormalization group equation ( |5^ in the sense of a formal power series in g iff it 
satisfies 

\^L-^-(v<j>,^^'^V^^\<j>,L) = -IC^'^\<P,L) (71) 

to every order r > 1. Notice that /C*-'"^(0, L) depends on V^^\(f),L) to lower orders 
1 < s < r — 1 only. Supplied with initial data, (^) can be integrated and yields a 
recursion relation for the coefficients in (|69|). 



4.1.1 First Order 

To first order (^TJ) is a homogeneos scaling equation 

It tells that the evolution is a pure scale transformation to first order. Consider the 
case of 0^-theory. There the first order is given by a {/(^-vertex 

V«(0,L)=A(L)i|d^x0(x)^ (73) 

The reader is invited to add a mass and a wave function term if he wishes. This 
0^-vertex is a solution to the homogeneous scaling equation ([72|) provided that the 
coupling flows according to 

A_(4-D)|a(L)=0. (74) 

The first order evolution is therefore A(L) = L^~^A(1). The conclusion is of course 
that the 0*^-vertex is relevant in dimensions D < 4, marginal in D = 4, and irrelevant 
in D > 4. In three dimensions the coupling flows proportional to the scale L. It 
diverges as L — * oo. To first order one already sees that global perturbation theory 
is unsuitable to perform this limit. 



4.1.2 Second Order 

The inhomogeneous term due to the renormalization group bracket of two (Z)"^- vertices 
is computed to 

/C(2)(0, L) = \{Lf J d^xd^y|^0(x)0(y) Cix - y) v{x - yf + 
i0(x)20(y)2 6C(x - y) v{x -y) + ^cj^^xf cj^^yf 20C(x - y)|. (75) 
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Two rather obvious remarks are in place here. First, the bracket is polynomial if 
all lower order vertices are polynomials. Second, the highest power of fields in the 
bracket is the sum of powers of the fused vertices minus two. It is instructive to enter 
the second order equation ([71|) with an ansatz containing precisely the interactions 
present in (75), 

V(2)(0,L) = XiLfJ d^xd^|/|^0(x)</.(|/) vPix-y,L) + 

vf {x -y,L) + vf (x - L) I . (76) 

It yields a solution to (|7TD iff the kernels in the ansatz obey the inhomogeneous 
differential equations 

- - Vg(x, L) = -/Cg(x) (77) 

with 

)C^^\x) = C{x)v{x)\ ICf\x) = 6C{x)v{x), /Cf (x) = 20C{x) (78) 

and 

af=3(D-2), aP=2{D-2), = D - 2. (79) 

Eq. (^) is easily integrated. Let us perform the integration in momentum space 
and thereby prepare the ground for the local perturbation expansion below. Fourier 
transformation turns ( [77| ) into 

- ?S + P^} WJiP, L) = -lcS!ip) (80) 

withg 

jcf\p) =C^v^v{p), }Cf\p) = 6C^v{p), K!i\p) = 20C(p), (81) 

and 

af^ = 2D-Q, =D-4, 5^ ^ = -2. (82) 

■^Thc convolution of two momentum space functions is here defined as F G{p) — 
{27r)-D Jd^qF{p-q) G{q). 
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Notice that our second order ansatz has a prefactor A(L)^. Therefore the scahng 
dimensions here are those of the corresponding monomials relative to the that of 
two (/(^-vertices. Another way to write the momentum space equation ( pOD is 

{^~^^" V£^(^P> L)] = -}C?r!iLp). (83) 
In this form it is immediately integrated to 

vSfe L) = L«' V<? (|, l) - "^t-^^KiSitp). (84) 

From this expression we can learn that the evolution tends to forget the irrelevant 
part of the initial potential as L >> 1, whereas the relevant part of the initial 
potential is enhanced. We will not renormalize the second order flow at this instant. 
But let us remark that due to the damping of irrelevant initial data an ultraviolet 
limit depends only on non-irrelevant interactions in the bare potential. Furthermore 
we see that the the integral in (|8^) contains divergent terms in the non-irrelevant 
case as L — > cxD. Both powerlike and logarithmic singularities appear. 

A general solution of the second order equation can be composed of this particular 
solution and any solution of the homogeneous scaling equation 

- 5g + V(J(p, L) = 0. (85) 

In the sequel we will adopt the following point of view regarding this freedom. The 
first order interaction enforces recursively through the bilinear bracket a certain set 
of vertices at higher orders. The homogeneous equation (|85|) allows us to introduce 
further vertices by hand into the iteration at higher orders, which are not present in a 
minimal scheme. We will restrict our attention to the solution to the renormalization 
group equation with a minimal set of vertices. This solution is determined by the 
first order potential. We can think of the first order as the germ of the theory. 
Introduction of other vertices at higher orders might however offer an interesting 
way to mix models. 

4.1.3 Higher Orders 

If the first order is polynomial, so are all higher orders in the minimal scheme. In 
the case of 0^-theory the general form of higher order vertices is 

V(^) (0, L) = X{LY J2 TTTTr I ^""3:1 . . . d^X2n 0(X1) . . . <P{x2n) (Xl, . . . , X2„, L). 

n=l l^^J- ■) (86) 
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This general form iterates to every order of perturbation theory. Other interactions 
could however be introduced by hand. The highest connected vertex built from r 
first order 0^- vertices has 2(r + 1) fields. The renormalization group equation for 
the kernels in (E&j reads 



^7) 



with real space scaling dimensions (T2^^ = m(2 + D) — r(4 — D). The equivalent 
equation in momentum space is 

2n-l 1 _ 

- ^2n + Pitt r "^^niVl, ■ ■ ■ ,P2n-l, L) = -/C2(pi, . . . ,P2n-l) 

dpi I (88) 

with momentum space scaling dimension = D + n(2 — D) ~r{4 — D). A compact 
way of writing (pSD then is 




which is then integrated to 



(89) 



L'SvS (^.■■- .^.l) + II, Jt-'''■^f:>P^,■■■,tp.,..^). (90) 

Thereby we have put perturbation theory for the renormalization group evolution 
into the form of a recursion relation. A recursion step consists of all contractions of 
previous vertices with one another to a given total order plus integration of (90). 

As it stands this perturbation expansion develops singular coefficients when the 
scale L is taken to infinity. These singularities can be removed by renormalization. 
An elegant way to renormalize the series goes as follows. In the recursion formula 
(90) the integral is performed from the ultraviolet end to the infrared end of the 
theory. Polchinski ||P84|| integrates the non-irrelevant degrees of freedom precisely 



the other way using a mixed boundary value problem. The relevant data is there 
prescribed on a lower scale than the irrelevant data. We mention that this approach 
can be generalized to a version of Symanzik's improvement program ||Wi88|| . This 
renormalization technology is well developed by now. 

In this paper we choose another route. We will search (and find) a way to 
formulate the theory in terms of quantities which are independent of the scale L. 
These quantities are auto-renormalized. They are identical with their scaling limits. 
The idea is switch to another form of perturbation theory. 
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4.2 Local Perturbation Theory 

We speak of local perturbation theory when the expansion parameter g is taken to 
depend on the scale L. A particular form of local perturbation theory is the run- 
ning coupling expansion for renormalized trajectories proposed in ||Wi96|| . There the 
ambition is to find renormalization group flows whose scale dependence comes exclu- 
sively in form of a running coupling g{L). The expansion coefficients themselves are 
scale independent. We will develop such an expansion for the 0^-trajectory in the 
infinitesimal renormalization group setup. A forerunner with discrete renormaliza- 
tion group transformation is found in |[Wi96|| . We intend to solve the flow equation 
in terms of a power series 



V{<P,g{L)) = j:^-^V^'H<P)- (91) 

r=l 

The interaction V^'^\(f)) is assumed to be independent of L. We choose the 0^- 
coupling as expansion parameter. It is deflned through the condition 

V,{0,0,0,giL)) = g{L) (92) 

on the four point kernel at zero momentum. To be precise, we impose the pertur- 
bative constraint that (^) be zero for all orders larger than one. The 0^-coupling 
is by no means the only possible choice in this approach. It is however a natural 
candidate when dealing with 0^-theory. It will be used in the following. In order to 
solve (0) in terms of (^TJ) we also require another power series 

L^9iL)=/3igiL)) = ±^-^/3^^\ (93) 
dL ^ r! 

The outcome of this analysis is a curve in the space of effective interactions para- 
metrized by g. The /3-function (^) says how the renormalization group acts on 
interactions on this curve. The coefficients Z?'-^^ form a second set of unknowns, be- 
sides those in V^'^''{(f)). They have to come out of the theory. To be more precise, 
the /3-function substitutes those degrees of freedom which are removed by (p2D. The 
intent with (|9l|) is a renormalization group orbit which is not parametrized by L 
but rather by g. As we will see such a curve is indeed determined to all orders of 
perturbation theory once we supply an appropriate first order interaction. In the 
case of 0^-theory the appropriate first order interaction is 

V(^)(0) = ^/d^^0(^)', (94) 

a 0^- vertex. The corresponding renormalization group orbit is called the 0*^-trajectory. 
It is the object of principle interest in massless ^'^-theory. Potentials on the 0^- 
trajectory are said to scale. The first order interaction (|9l) turns out to require a 
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slight modification in four dimensions. We need to add also a wave function term. 
We will do this below. 



4.2.1 Scaling Equations 

We insert the expansion (^) into the renormalization group equation (Q) and 
deduce a system of equations for the unknowns therefrom. This system is called the 
set of scaling equations. We have 

Therefore, our local perturbation expansion is a solution to (|54D in the sense of a 
formal power series iff the unknowns obey 

r / \ ^ X \ r — l / \ 



(96) 



to all orders r > 1. The most remarkable property of (^) is again independence 



of L. Another way to think of (|96D is that we are looking for a fixed point of the 
transformation composed of a renormalization group step and a transformation of 
the four point coupling. We organize (0) into a recursion relation. 

4.2.2 First Order 

Scaling requires of V^^\(p) to be a scaling field of the trivial fixed point. To first 
order (|96|) reads 



|/?«-(l)0,A^|v«(0) = O. (97) 

It requires of V^^-* (0) to be an eigenvector of the generator of dilatations, and of P^^^ 
to be its eigenvalue. The ^'^-vertex is indeed an eigenvector. Its eigenvalue has 
the familiar value 

=4-D. (98) 

Here three remarks are in place. First, if we want to add further terms to the 
first order (|94D then (^) requires of all other terms to be eigenvectors with this 
same eigenvalue (|9^). In four dimensions the eigenvalue ( pH] ) is zero. There one 
has indeed another marginal eigenvector, the wave function term. Second, we could 
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have considered any other scahng field of the trivial fixed point as our first order 
starting point, for instance a 0^-vertex in three dimensions. All of the scaling fields 
generate interesting trajectories, for which this theory equally well applies. Third, 
the eigenvalue does not need to be larger or equal to zero. We can as well start 
with an irrelevant scaling field. In fact the 0^-vertex is marginally irrelvant in four 
dimensions. 



4.2.3 Second Order 

The first order equation (|97|) is special in that it is homogeneous. All higher order 
scaling equations are inhomogeneous. The second order equation is given by 

|2/3(i) - I V^'\<P) = -/3(^)V(1)(0) - 2 (V(i), V(i)) . (99) 

Notice that unlike in the global expansion there is as well an inhomogeneous term 



due to the fiow of the coupling parameter. Minus the right hand side of (|99|) is 
explicitely computed to 



2 r ^x)(j){y)C{x - y)v{x - yf + 



j^(t>{xfcP{yfQC{x - y)v{x - y) + ^^<P{xf <P{yf2QC {x - y)\. (100) 



The general solution of ( p9D consists of special solution plus any solution of the 
homogeneous equation. In the global expansion special solutions were singled out 
by boundary conditions. Here we will do something different. We will only admit 
solutions which are given by finite kernels in momentum space. Furthermore, we 
will restrict our attention to local interactions in the minimal scheme. It then turns 
out that there exists only one such solution. This special solution contains precisely 
those vertices, which are enforced by the first order interaction. Therefore, we make 
the following ansatz 

= E tAt / d^'xi . . . dPx2n </)(a;i) . . . (t){x2n) V£^(xi, . . . , X2n). 

n=A^^)--' (101) 

From (100) it is clear that we can only succeed in terms of distributional kernels. 
We therefore Fourier transform (PP|) and thereby obtain 



2n-l Q ~l 



2/5(1) -{D- n{D - 2)) + ^ p^-— V£)(pi, . . . ,P2n-l) = -4n^(pi, . . . ,P2n-l) 

m=l '^Pm) (102) 
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with inhomogeneous terms given by 

K:P{p) = Ci.vi.v{p), (103) 
jcf\puP2,P3) = f3^^^ + 6C^v{p,+p2), (104) 
4'^(pi,... ,P5) = 20C(pi+p2+P3). (105) 

The right hand sides are here understood to be symmetrized in the momenta. We do 
not write this symmetrization exphcitely in order to simphfy the notation. Here we 
recognize again the the second order scahng dimension (T2n — 

which was introduced above. The general method to solve the differential equation 
( |102| ) is explained below. We immediately apply it to this second order equation. 



Consider first the six point interaction. Since cxg = —2, it is irrelevant (in any 
dimension), we can immediately integrate the six point interaction toQ 

(pi, . . . ,p5) = - /' ^t'JC^\tp„ . . . ,tp,). (106) 
Jo t 

Notice the similarity with the analogous expression in the global expansion. For- 
mally it can be obtained by taking L to infinity in (90). We then treat the four 
point interaction. It is computed in two parts. The first part concerns its value at 
zero momentum. By definition of our expansion parameter we have to satisfy 

= /?(2) + 6(7^^(0). (107) 

Using the (/)^-coupling to organize the expansion, we imposed that V4^'*(0, 0, 0), and 
all higher orders of the four point kernel at zero momentum, be zero. Eq. (|107|) 
determines the value of In four dimensions its value is computed to 

As a consequence the 0^-vertex is indeed marginally irrelevant in four dimensions. 
Thereby the renormalization group flow on the 0^-trajectory is not asymptotically 
free in the ultraviolet. Nevertheless the 0^-trajectory is a well defined object at weak 
coupling. The irrelevant remainder of the four point vertex is then integrated to 

Vf (pi,P2,P3) = -6 jt^'''{C^v{tpi+tp2)-C^v{0)}. (109) 



^It is instructive to perform this integral. The result is Vg^^(pi,... ,^5) = 10(pi + P2 + 
Pa)"^ {exp {—{pi + P2 + Ps)^) ~ l}- This expression is regular at zero momentum and a bounded 
function of the momenta. It is of the form of a cutoff propagator. 
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Here it becomes transparent how the flow of the couphng parameter saves us from 
a logarithmic singularity in four dimensions. Notice also that the subtraction of 
the zero momentum piece is unnecessary in three dimensions. There this vertex is 
already irrelevant. The two point kernel finally requires most attention and poses 
even an obstacle at this stage. Due to Euclidean invariance we have 

W\p)=A{p'), K,^?{p)=2B{p'). (110) 

It is convenient to trade for a new variable u. The scaling dimension of the two 
point kernel is a^'^ = 2D — 6, and is two in four dimensions. The scaling equation 
for the two point kernel then becomes 

- l] A{u) = -B{u). (Ill) 




Its solution requires a second order Taylor expansion with remainder term for the 
function 



„2 ,1 



A{u) = A{0) + u A'(0) + — / dt(l - t)A"{tu). (112) 

2 JO 

The zero momentum value is directly determined. Its value is 



A{Q) = 5(0) = ]^C^v^ m = ^ (2 log(2) - log(3)) , D = A. 



(113) 



It can be interpreted as a second order mass parameter. The Taylor remainder in 
( |112| ) is computed from the second derivative of ( |111| ). We have 

\^u^ + l^A'\u) = -B'\u). (114) 

Two ^-derivatives have changed the powercounting by four units. As a consequence 
the second u-derivative is irrelevant and is therefore integrated to 

A"{u) = - f' dt B"{tu) (115) 



in complete analogy to the case of the six point kernel.^ The grain of salt is the first 
M-derivative, a wave function term. Its scaling equation is 

u^A'(u) = -B'iu). (116) 

mi 



^ A parameter representation for the inhomogeneous mass term is /Cj^^ (p) — C -k v -k v{p) — 
2(47r)^-^ dai da2{aia2 + ai+ a-})^^!'^ exp {~a\a2'p^ j (a\a2 + ai + ct-ii) ■ Notice that it is 
regular function of the variable in four dimensions. 
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Since it is marginal, it does not determine the first w-derivative at zero momentum. 
Furthermore, it requires 



B'{0) = 0. (117) 

Eq. ( |117D can be checked to be false in four dimensions. The problem is a wave 
function term which is generated dynamically to second order. As it stands the 
theory is inconsistent. This problem is solved by introduction of a wave function 
term to first order. Thus we extend (|9^) to 

V(i)(0) = ^ / d^x(j){x)^ + ^ / d^xd^(f){x)d^(j){x). (118) 



4! 7 ' 2 

A first order wave function term is consistent with (^) in four dimensions. Both 
terms in ( p.l8| ) are marginal. The first order scaling equation leaves the new param- 
eter C^^) undetermined. We use this freedom to satisfy the second order constraint 
( |117| ). It is a general feature of this approach that marginal parameters are deter- 
mined one order later than the other ones to a given order. Nevertheless we have a 
recursive perturbation theory. It is now consistent to all orders. Having accompa- 
nied the 0^- vertex with a first order wave function term, we find a few more effective 
interactions in the analogue of (100). With (0) replaced by ( |11^ ), it becomes 



{x)(t){y)C{x - y)v{x - yf + 



^0(x)V(z/)'6C(x - y)v{x -y) + ^^<P{xf (t){yf2QC {x - 1/)| + 
/^^'^C^'^ Wxd,<P{x)d^<P{x)+ /d^a;d^y|^^0(a;)0(y)2(-A,)(-A,)C(x-y) + 



^^<P{xf2{-Ay)C{y)v{y) + ^^^{x)<P{yf%{-l\y)C{y). (119) 

It follows that we have to replace eq. (105) by 

^f(p) = C*ty^tJ(p) + C^'^a + /3^'^C^V + 2(C^'^)'(p')'(5(p), (120) 
(Pi,P2,P3) = + QC^v{p^+p2) + 8C(') plC{pi), (121) 

(Pl, . . . ,P5) = 20^(^1+^2+^3). (122) 

The constant a is here given by a the convergent one loop integral 

d^p 



2 f d''y{-Ay)C{y)v{y) = 4 / ^exp(-2p2) = (47^)-^ D = 4. 
J J [ZTT) ^^23) 
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We then proceed exactly as above. The new terms in (122) do not disturb the 
recursion. The six point kernel is not affected at all. Neither is the four point kernel 
at zero momentum. Therefore the coefficent /J'-^-* is independent of the new constant. 
The quartic remainder gets an extra contribution and requires first knowledge of 
Spelled out explicitely, the eq. ( p.l7| ) becomes 



= P^'\'^'^ + ^^C^v^v{p) 



(124) 

p=0 



It determines the first order wave function parameter. The second order work is 
therefore best organized as follows. We first compute Z?*^^-* from (|107|) , and second 



(^(1) from ( |124| ). Knowing of the marginal data we then compute all kernels in terms 
of their Taylor expansions and remainders. In the normal ordered formulation the 
order in which the kernels are computed is of no importance. The value of the first 
order wave function paramter comes out as 

C^'^ = 7\, - (125) 
^ 18(47r)2 ^ ' 

The second order effective mass parameter is then changed from (|113|) to 



^(0)=(i)l(2M2)-log(3)-l). (126) 

The changed remainder terms of the quadratic and the quartic kernel follow immedi- 
ately from (122). This completes the calculation of the second order. In summary we 
have extracted the non-irrelevant part by Taylor expansion. The Taylor coefficients 
have been determined directly by evaluation at zero momentum. The irrelevant 
part, including the Taylor remainders, have been obtained by integration of the 
corresponding scaling equations. 



4.2.4 Higher Orders 

The second order scheme generalizes to third and higher orders. Assume that we 
have computed V^'^\(f)) and P^^^ to all orders 1 < s < r — 1 except for (^('~~^). Then 
we first compute P^'^\ thereafter C^^~^\ and then the V^'^\(f)) except for The 
non-irrelevant degrees of freedom are again separated out by Taylor expansion in 
momentum space. A convenient notation for this Taylor expansion goes as follows. 
Introduce scaling fields 

02,o(0) = 11^''^ H^f^ ^2,2(0) = ^ / d^'x d^cl>{x)d,^{x), (127) 
(^4,o(0) = ^/d^x0(x)^ (128) 
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with scaling dimensions 



(V(j), j-] 04,0(0) = (4 - D)02,o{(p). 



(129) 
(130) 



In four dimensions, the non-irrelevant part of the effective potential can be written 
as a sum 

VS(0) = f^^'^O^M + C^^-^O^M) + A^^^04,o(0). (131) 
The general form of the effective potential to order r in the minimal scheme is 

„=i(2n)!^ (132) 
The non-irrelevant coupling constants are, in terms of these kernels, given by 



2D 



d^X2 (X1-X2)' V('-)(Xi,X2) 



9 



9(p2 



= J d^X2 d^X3 d^X4 VP(X1,X2,X3,X4) = (0,0,0). 



5 




(133) 


yvr^(p) 




(134) 


(0,0,0). 


(135) 



A compact notation for their extraction from the polynomial (|132|) is the formal 
pairing 

/x^'-) = (02,0, V^'-)) , C^'^^ = {02,2, V^'-)) , A^'^) = (04,0, VW) . (136) 

Its meaning is, as was said before, nothing but Taylor expansion in momentum 
space. The scaling equation to order r can be put into the form 



6 



(137) 



The right hand side of (|137|) still depends on /^^^^^ and C*-*^ It is given by 

/CM(0) = ^ (''\ /5Wv(^-^+i)(0) +'^' ("^ (y^'\V^'-'^) . (138) 



s=l 
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Projecting both sides of ( |137| ) to 04^0 (0) "we find 

i)/?(i)am = - J2 n /5('^)A(^-^+^) - J2 i"^] (^4,0, v^'--^) 



r 



s=2 \ / s=l \ / 



(139) 

Since A^^-* = by construction, eq. ( |139| ) determines the coefficient The right 
hand side of ( [1^ does not depend on C^''"^^ since (C4,o, (v(^)(0), 02,2)) = 0. Thus 
the order r — 1 wave function term does not contribute through the renormahzation 
group bracket to an effective ^'^-vertex. Thereafter projecting both sides of ( |137| ) to 
02,2(0) we find 

In four dimensions the left hand side is zero since 

= 4 - D = 0. But then 

( |14U| ) determines the value of C}~~^~^^ . Notice that Z^^^-* is not zero. Notice fur- 
ther that the renormahzation group bracket does not depend on C'''^"^'' because 
(^^2,25 02,2)) = 0. The reason is that the covariance C is regular at zero 
momentum. The rest of the work is immediately put to order. The effective mass 
parameter to order r follows from 



(r/3(^) - 2) /iM = - E ( ' ) /^^^V^^-^^^^ - E ( ' ) ((^2,0, (V(^), V(^' 



s=2 



(141) 

Its computation requires the extraction of the effective mass term in the renormal- 
ization group bracket. The computation of V^g](0) except for C^^ is complete. The 
irrelevant part is directly integrated following exactly the scheme of the second order 
calculation. Eq. ( |138| ) can be expanded into a polynomial 

^^'^0) = E / d^'xi . . . d^X2„ 0(Xi) . . . 0(X2„) /C(2'-)(xi, . . . , X2„). 

n=l l^^J- (142) 

The scaling equations for the Fourier transformed kernels are given by 

\ -^2n + E P^TS r ^^n{Pl, ■ ■ ■ ,P2n-l) = -^tl{PU ■ ■ ■ ,P2n-l)- 

For n > 3, the scaling dimension is a2n = D + n{2 — D) — r(4 — D) = A — 2n < 0, 
D = 4. Those kernels are therefore all irrelevant. The scaling equation is in this 
case integrated to 

VSiPU . . . ,P2n-l) = - / -t-"^- ^tnitPl. ■ ■ ■ ,tp2n-l). (144) 

Jo t 
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The integral converges because of the negative power counting. The zero momentum 
part of the four point kernel has already successfully been transferee! to the j3- 
function. Its remainder is irrelevant and integrated to 

vt\pi,P2,P^) = - I' J {ict\tp„tp„tps) - )C^\0,0,0)} . (145) 

The integral converges due to the subtraction at zero momentum. Finally, the two 
point kernel is reconstructed with the help of 

Wip) = A{p'), )C^\p) = 2B{p'), (146) 

and 

A{u) = /i^'^) + C^'^M + — ds(l - s)A"{su) (147) 

2 JO 

through 

A"{u) = ~ [\tB"{tu). (148) 

JO 

The scheme is now complete. We have a manifestly finite recursive local perturbation 
theory. 



5 Regularity 

The integration of renormalization group differential equations generally requires 
initial data. In the local expansion, where we no more have evolution equations, we 
substitute the initial data by a requirement that the solutions be finite and regular. 



5.1 Renormalization Group PDEs 

In the perturbation expansion the recursion to each higher order consists of solving 
a set of renormalization group PDEs of the general form 

p^-a^F{p)=G{p). (149) 

Here a G Z, and G{p) is a given function of p G M^. In this section we will have a 
look at the regular solutions of ( [149| ) . 



31 



5.1.1 Irrelevant Case 



The irrelevant case is defined by a < 0. Let us assume that G{p) is a continuous 
function on M^. Eq. ( |149| ) is equivalent to 

L-^ {L-^'FiLp)} = L-"G{Lp). (150) 



A special solution to ( |149|) is then given by 

F{p)= r^L--G(Lp). (151) 

Let us require the solution to be a continuous differentiable function on M^. Then 
we notice the following facts: 

I) There exists a unique solution to II) It is given by ( \151\ ). 



First, suppose that we have two different solutions Fi{p) and F2{p) of ( p.49| ) which 
are both continuous differentiable. Their difference satisfies 

^ -a]{F,{p)-F2{p))=0. (152) 




The only solution to this equation, which is regular at the origin in the case ex < 0, 
is zero. Second, ( |151j ) is continuous differentiable and is a solution to (|149|) . 



Thus we can substitute initial or boundary data by regularity to obtain a unique 
solution. Notice that its value at zero is — crF(O) = G{0), as can be seen from both 
and (imp. 



5.1.2 Relevant Case 

The relevant case is defined by o" > 0. In this case we cannot use (|151|) because the 



integral diverges, unless G{p) provides a sufficiently high power of p. The trick is 
to perform a Taylor expansion with remainder term to high enough order. Let us 
assume that G{p) is (a + l)-times continuous differentiable. 

Let a = («!,... ,aAr) G be an integer valued multi-index. Define |a| = 
«! + • ■ ■ + Oiv, a\ = ail . . . Oat!, and = Pi^ . . -p^^ ■ Then F{p) is a solution to 
( |149| ) iff its derivatives with a < |q;| satisfy 

Thus each momentum derivative reduces the power counting parameter by one unit. 
Solutions to ( |153| ) which are regular at the origin satisfy 

-(^-H)^(0) = ^(0)- (154) 
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Their Taylor coefficients are all determined except for those with a — \a\ = 0, the 
marginal ones. Let us assume that 

^(0)=0, a-\a\=0. (155) 
Then ( |149[ ) can be solved in terms of a Taylor expansion of order a with remainder, 



a a\a\ p a f)\a\ p 

nrt- E ^-^(0)+ E d'd-T-^Crt. (156) 

|cj|<cr \a\=a+l 

The derivatives of order a + 1 have negative power-counting and are integrated as 
above. Let us require the solution to be [o + l)-times continuous differentiable. 
Then we have: 



I) There exists a set of solutions to (\14^ which can be parametrized by its Taylor 



coefficients with a = \a\. II) The relevant Taylor coefficients with a < \a\ are 
uniquely determined by l \154 )- HI) The Taylor remainder is reconstructed from the 



{a + l)th derivatives. It is unique and follows from (\151\) 



5.2 Large Momentum Bound 

We prove a large momentum bound for the solution to the renormalization group 
PDE under the assumption of a large momentum bound on the inhomogeneous side. 
We choose an Loo,£-norm for some e > 0. It is rather wasteful but suffices to prove 
finiteness of the bilinear renormalization group bracket. 

5.2.1 Irrelevant Case 

Let a < 0. Suppose that the function G{j>) in (|149|) has a finite Loo,e-norm 

||G|U,e = sup [\G{p)\e-'\p\] < oo. (157) 



Then the solution ( |151D inherits an Loo,e-bound. From 

"1 dL . , /-i dL 



|F(p)|e-^|p| < / — L-"|G(Lp)|e-^IPl < / — L-'e-^^-^WPl ||G|| 
Jo L Jo L 



oo,t 

(158) 



it follows that 



< — IIGIloo,.. (159) 
—a 



Eq. (|159|) shows that the irrelevant solution to the renormalization group PDE is 
not only finite but also decreases in the Loo,e-norm. 
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5.2.2 Marginal Case 

Let a — 0. In this case we assemble F{p) using a first order Taylor formula. Suppose 
then that we have Loo,e-bounds on the first derivatives 



\G 



//||oo,e 



sup 



d 



G{p) 



(160) 



If F{p) is marginal, then its first derivatives are irrelevant with scaling dimension 
minus one. It follows that 



IF II < lie II 

I""- fj,\\cx),e — ||^/i||oo,e- 



Therefrom it follows that 



\F{p)\e-'\^\ < \F{0)\e-^\P^ + J2\P,\ [' dt\F,{tp)\e 
<\Fm + J2\p,\ /'dte-(^-*)^l^l||F^| 



(161) 



(162) 



The result is an Lr^ ^-bound 



\F\U<\Fm + -^\\F,\ 



(163) 



This estimate is not uniform in e. It works for e arbitrary small, but the bound 
grows with an inverse power of e. The large momentum growth is a consequence of 
the split in derivatives and Taylor remainder. 

5.2.3 Relevant Case 

Let (J > 0. This case requires a generahzation of the bound in the marginal case. The 
Taylor expansion is pushed to order a + 1. Then the derivatives become irrelevant. 
We assume ^-estimates on all derivatives 



IG'allco.e — sup 



G{p) 



-e|p| 



< OO 



(164) 



of order \a\ = a + 1. Since they are irrelevant with scaling dimension minus one it 
follows that the corresponding derivatives of F(p) obey 



IF II < IIG II 

K alloc, e _ ||^a||oo,e) 



(165) 
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and are also Loo,e-bounded. Prom the Taylor formula it then follows that 

|a|<<T |Q:|=fT+l " 

|a|<(T |a|=a+l (166) 



and thus 



ll^lloo,e< E ^^^,Hl^a(0)|+ E ^^.,a+l|l^a||oo,o (167) 
|a|<(T ■ |a|=(T+l 

with constants 

= sup {|p|l"'e-l^l} , = 3^^. (168) 

Thus we again have an L(x,,e-bound on the function F{j)). This completes the large 
momentum bound on F(p). Exactly the same strategy applies to the derivatives 
of Fiji) as well. The irrelevant derivatives inherit immediately large momentum 
bounds. The relevant derivatives require Taylor expansions. We omit to spell out 
explicitely the necessary bounds on the derivatives of G(p). 



5.3 Iteration and Regularity 

The iterative scheme determines order by order C*"*^^"*) ^"^^ the irrelevant 
remainders Vf^,^2n(Pi) • • • )P2n-i)- It is finite to all orders of perturbation theory 
because of the following iteration of regularity. Suppose that we have shown the 
following to all orders s < r — 1: 

I) f3^^\ C^^~^\ and n^^^ are finite numbers. II) Virr,2n(Pi) • • • iP2n-i) is a smooth 
function on IR x ••• x M for all 1 < n < s -\- 1, symmetric in the momenta, and 
0{D) -invariant. Ill) || Virr,2n,Qlloo,e is finite for a// e > 0, 1 < n < s + 1, and \a\ > 0. 
Here a is a multi-index which labels momentum derivatives. 

Then the same statements hold at order s = r. Since they are trivially fulfilled to 
order one they iterate to all orders of perturbation theory. 

To prove the iteration of regularity we once more inspect each step of the iterative 
scheme. First, the irrelevant remainders ^lrr,2n(Pi) • • • )P2n-i) are smooth functions 
on IR X • • • X M, symmetric under permutations and 0(D) -invariant. They and 
all their momentum derivatives satisfy Loo,e-bounds. They are composed of two 
contributions. The first immediately inherits a bound from the induction hypothesis. 
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The second is a sum of renormalization group brackets of lower orders. Therefore, 
they consist of multiple convolutions with propagators. The integrals converge, are 
smooth functions of the external momenta, and satisfy Loo,e-bounds. Second, we 
have linear equations for the coefficients P^^\ C^^~^\ and fi^^"^ with finite coefficients. 
Third, the integration of the inhomogeneous renormalization group PDEs, yields 
solutions with the desired properties. 



6 Conclusions 



The aim of perturbative renormalization theory is to derive power series expansions 
for Green's functions which are free of divergencies. The BPHZ theorem states that 
this can be accomplished by writing the Green's functions in terms of renormalized 
parameters. An elegant proof of the BPHZ theorem was given by Callan [ |G76| ]. A 
polished version of which is due to Lesniewski [[L83|| . Their method is similar to ours 



in that it is based on renormalization group equations for the renormalized Green's 
functions, the Callan-Symanzik equations. In some sense (|I|) is a Wilson-analogue 
of the Callan-Symanzik equations. The method proposed here is different in that it 
does not resort to any kind of graphical analysis, not to analysis of sub-graphs, and 
not to skeleton expansions. 

A new generation of proofs of the BPHZ theorem was initiated with the work 
of Polchinski ||P84]| . His proof has been simplified further by Keller, Kopper, and 
Salmhofer | |KKS9(J[ | . Their approach is similar to the method advocated here in 
that it is based on Wilson's exact renormalization differential equation. The details 
are however quite different. The main difference is that Polchinski begins with a 
cutoff theory. He then shows how the cutoff can be removed in a way such that the 
effective interaction remains finite. Our method directly addresses the limit theory 
without cutoffs, expressed in terms of a renormalization group transformation with 
cutoffs. In some sense we are here simultaneously changing Polchinski's renormal- 
ization conditions and integrating an amount of fluctuations. Unlike Polchinski and 
followers we use a renormalization group differential equation with dilatation term. 
A way to think of (|I|) is as a renormalization group flxed point of a system which 
has been enhanced by one degree of freedom, the running coupling. This flxed point 
problem can only be formulated with rescaling and with dilatation term. 

Another renormalization group approach to renormalized perturbation theory 
comes from Gallavotti |p85| , |GN85|| and collaborators. Pedagocial accounts of tree 
expansions can be found in |[BG95| , |I''HRW88| . There the result of renormalization 



is expressed in terms of a renormalized tree expansion. The program of |[Wi96|| with 
an iterated transformation with flxed L is related to the tree expansion. Both are 
built upon a cumulant expansion for the effective interaction. The renormalization 
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procedure is however quite different. Like Polchinski, Gallavotti starts from a cut- 
off tlieory. It is organized in terms of trees, wliicli describe tlie sub-structure of 
divergencies in Feynman diagrams. Tlie divergencies are transformed into a flow 
of tlie non-irrelevant couplings. This part is similar to ours. The basic difference 
with Gallavotti is that we do not organize our expansion in terms of trees. A hybrid 
approach between Polchinski and Gallavotti is due to Hurd | li89 |. 

An important question is wether this construction of renormalized trajectories 
extends beyond perturbation theory.^ Another important question is wether it ex- 
tends to renormalized trajectories at non-trivial fixed points. We hope to return 
with answers to these questions in the future. 
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